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1.0  INTRODUCTION 

The  evolution  of  nonintrusive  flow  diagnostic  methods  in  the  last  few  years  has  created 
powerful  tools  for  gathering  meaningful  flow-field  data  in  extremely  complex  flow  fields.  But, 
as  with  any  measurement  technique,  care  must  be  taken  in  defining  the  limitations  and  critical 
requirements  for  these  techniques.  This  is  especially  true  in  the  case  of  these  nonintrusive 
techniques  since  there  are  often  no  means  to  verify  measurements  in  the  complex  flow  fields 
in  which  the  techniques  are  being  applied. 

One  of  the  most  promising  of  these  techniques  is  laser  velocimetry  (LV).  The  LV  measures 
the  velocity  of  small  particles  entrained  in  the  flow  field  of  interest.  The  fluid  velocity  and 
turbulent  velocity  correlations  are  then  inferred  from  statistical  analysis  of  these  particle  velocity 
measurements  under  the  assumption  that  the  particles  exactly  reproduce  the  motion  of  the 
fluid.  This  assumption  needs  to  be  carefully  evaluated  since  it  is  fundamental  to  the  application 
of  the  laser  velocimeter. 

Early  development  of  the  velocimeter  seemed  to  have  centered  on  the  electronics  and  optics 
of  the  instrument.  Particle  dynamics  was  generally  assumed  to  be  a  minor  problem  compared 
with  the  task  of  acquiring  any  data  at  all.  Simple  low-speed  flow  fields  were  used  to  verify 
the  LV  principle,  and  calculations  using  simple  particle  drag  laws  and  simplified  particle 
equations  of  motion  were  used  to  dispel  concerns  about  particle  dynamics  effects.  As  the 
LV  was  graduated  to  more  complex  and  higher  speed  flows,  it  became  obvious  that  the  assumed 
small  particles  were  not  following  the  fluid.  This  can  be  seen  in  the  exit  of  a  supersonic  jet 
as  described  in  Ref.  1.  Results  such  as  these  have  led  to  a  more  thorough  evaluation  of  particle 
dynamics  and  particle  generation  in  the  LV  community. 

Most  of  the  current  particle  dynamics  studies  are  directed  to  the  problem  of  determining 
an  accurate  mean  flow  velocity  measurement.  Very  little  has  been  or  is  being  done  to  investigate 
particle  response  to  the  velocity  fluctuations  of  a  turbulent  flow  field.  Hjelmfelt  and  Mockros 
(Ref.  2)  performed  a  fairly  simple  analysis  of  particle  response  to  a  sinusoidally  oscillating 
flow  field.  The  solution  is  only  of  academic  interest  since  turbulence  consists  of  finite  random 
eddies  rather  than  continuously  varying  velocity  oscillations.  The  solution  is  based  on  the 
questionable  assumption  that  the  solid  particle  always  follows  the  same  fluid  particles;  that 
is,  there  are  no  trajectory  crossing  effects.  Sheun  et  al.  in  Ref.  3  modeled  the  turbulent  field 
more  accurately  by  using  a  Monte  Carlo  approach  to  reconstruct  the  randomness  of  turbulence 
and  included  trajectory  crossing  effects.  This  technique  is  adequate  for  large,  heavy  particles 
but  lacks  the  detail  in  the  turbulence  model  and  in  the  particle  equation  of  motion  to  adequately 
predict  the  motion  of  the  small  particles  of  interest  in  velocimetry  applications.  Ormancey 
and  Martinon  (Ref.  4)  followed  a  similar  approach  to  Shuen  but  included  a  more  complete 
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model  of  a  turbulent  eddy  and  an  improved  equation  of  motion.  This  method  fails  to  predict 
the  motion  of  small,  light  particles  in  a  turbulent  flow,  which  possibly  indicates  some  problems 
exist  with  his  turbulence  model.  Since  small,  light  particles  that  nearly  reproduce  the  motion 
of  the  turbulent  fluid  are  of  principal  interest  in  LV  applications,  improvements  to  these  models 
are  necessary  before  they  are  useful  in  evaluating  particle  dynamics  limitation  for  LV 
measurements. 

The  objective  of  this  effort  is  to  develop  the  models  necessary  to  determine  particle 
requirements  for  obtaining  turbulence  measurements  with  the  LV  and  to  apply  these  models 
to  some  simple  flow  fields  to  evaluate  particle  response.  The  particle  equation  of  motion  is 
carefully  evaluated  in  Section  2.0,  and  a  computer  code  is  developed  to  solve  this  equation. 
Using  this  code,  the  broadening  of  the  velocity  distribution  attributable  to  a  particle  size 
distribution  in  the  presence  of  a  velocity  gradient  is  discussed  in  Section  3.0.  A  Monte  Carlo 
turbulence  model  is  developed  and  coupled  with  the  particle  equation-of-motion  solver  from 
Section  2.0  in  Section  4.0.  This  code  is  applied  to  two  simple  turbulent  flow  fields  in  Section 
5.0  to  evaluate  the  particle  size  requirements  for  turbulence  measurements  in  these  flows.  Simple 
“rule  of  thumb”  techniques  for  estimating  particle  size  requirements  for  turbulence 
measurements  with  the  LV  are  outlined  in  Section  6.0.  Conclusions  and  recommendations 
are  presented  in  Section  7.0. 

2.0  PARTICLE  EQUATION  OF  MOTION 

There  are  a  number  of  forms  of  the  equation  of  motion  of  a  small  particle  in  a  nonuniform 
flow  field  in  the  literature.  All  these  equations  are  based  on  the  equation  of  motion  for  a 
small  sphere  accelerating  in  a  fluid  at  rest  developed  independently  by  Basset  (Ref.  5), 
Boussinesq  (Ref.  6),  and  Oseen  (Ref.  7).  Tchen  (Ref.  8)  generalized  the  solution  for  unsteady 
Stokes  flow  in  a  uniform  fluid.  Tchen  then  attempted,  through  heuristic  argument,  to  extend 
the  solution  to  nonuniform  flows.  It  is  assumed  in  all  these  equations  that  the  particle  is  small 
enough  to  not  disturb  the  flow  field  in  which  it  is  contained.  This  leads  to  the  requirement 
that  the  particle  loading  be  extremely  low. 

2,1  BASIC  EQUATION 

Tchen’s  equation  may  be  written  as 

dvj  Duj  1  d  , 

mp  — —  =  -  3'/r/t  d(v;  -  us)  +  mf  -  —  mf  —  (vj  -  Uj) 
dt  Dt  2  dt 
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where 


and 
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'J'TQfH  j 
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d  1 

a;  <v>  ~  u'> 

sjt  —  T 


dr 


-  a  (mp  -  mf) 


(1) 


Duj  3uj  3u, 

-  =  - — -  +  u ; - 

Dt  3t  dxj 


(2) 


du;  dUj  dllj 

-  -  -  +  V; - 

dt  dt  dxj 


(3) 


and  Vj  =  v,(t),  u,  =  u;(t).  The  terras  on  the  right-hand  side  are  the  Stokes  drag  term,  resultant 
pressure  term,  apparent  mass,  Basset  history  integral,  and  the  gravitational  force.  The  terms 
in  Eqs.  (2)  and  (3),  which  were  not  rigorously  derived,  are  a  source  of  controversy  with  this 
equation. 

Several  authors  have  attempted  to  improve  upon  Tchen’s  equation  of  motion  for  a  sphere 
in  a  nonuniform  Stokes  flow,  including  Corrsin  and  Lumley  (Ref.  9),  Buevich  (Ref.  10),  Soo 
(Ref.  11),  and  Gitterman  and  Steinberg  (Ref.  12).  In  a  more  recent  effort  by  Maxey  and  Riley 
(Ref.  13),  the  equation  of  motion  is  derived  using  matched  asymptotic  expansion  of  the  inner 
disturbance  flow  and  the  outer  undisturbed  flow.  Their  derivation  also  included  the  effects 
of  velocity  curvature  on  sphere  drag.  Maxey  and  Riley’s  equation  may  be  written  as 
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dvj 
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+  (nip  -  mf)  gj 


(4) 


where 
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dt 


5uj 
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(5) 


dui 

dt 


3uj 

ir 


X  =  Y  (t) 


(6) 


and  Vi  =  Vi(t),  Uj  =  [  Y  (t),t].  Y  (t)  is  the  particle  position  as  a  function  of  time  in  an 
inertial  reference  frame.  Crudely  evaluating  the  order-of-magnitude  of  the  term  in  Eq.  (1) 
or  Eq.  (4)  relative  to  the  left  hand  of  the  equation  we  have 
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For  the  small  and  relatively  heavy  particles  of  interest  in  velocimetry,  the  drag  term  clearly 
dominates  the  equation  of  motion.  As  the  density  ratio  decreases,  other  terms  in  the  equation 
will  become  more  important. 

2.2  BASSET  TERM 

The  Basset  term  is  frequently  omitted  from  solution  schemes  for  the  particle  equation 
of  motion,  because  it  is  not  easy  to  deal  with  numerically.  Order-of-magnitude  analyses  similar 
to  those  in  the  previous  section  are  usually  used  to  justify  the  omission  of  this  term.  These 
analyses  fail  to  consider  the  magnitude  of  the  integral  term  and  should  be  viewed  with  caution 
if  the  magnitude  of  the  integral  is  not  known.  The  Basset  term  is  essentially  the  time  integral 
of  the  relative  acceleration  of  the  particle  weighted  by  the  square  root  of  the  time  difference 
between  the  present  time  and  the  time  at  which  the  acceleration  occurred.  It  is  significant 
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when  the  particle  relative  acceleration  is  large  compared  to  the  particle  relative  velocity.  It 
also  becomes  more  significant  as  the  particle  density  approaches,  or  is  less  than,  the  fluid 
density.  Tchen’s  solution  of  the  particle  equation  of  motion,  though  not  entirely  accurate, 
offers  the  opportunity  to  evaluate  when  the  Basset  term  becomes  important.  TWo  simple 
examples  using  Stokes  drag  law  follow. 

The  solution  to  Tchen’s  equation  of  motion,  including  all  terms,  for  a  step  change  in  fluid 
velocity  from  u  =  b(,  to  u  =  a0  is 


v 


e-kt  (b0  -  ac)  (cos  cot  +  —  sin  cot)  +  a0 

co 


(7) 


where 


-  k2  =  a2  |3  p  + 

> 

902  I 

4  | 

/ 

(8) 

k  =  a(l  -  3/2  j 8) 

(9) 

12  V  a 
a  -  - p 

(10) 

d2 

=  3(2  qp/q{  +  l)-> 

(11) 

If  the  Basset  term  is  ignored,  the  solution  becomes 

V  =  (bQ  -  aoje-"1  +  a0  (12) 

Solutions  of  these  equations  for  density  ratios  S  of  1000,  5,  1.75,  and  1.0  are  shown  in  Fig. 
1.  Fig.  la  shows  that  the  Basset  term  is  not  negligible  for  S  =  1000,  which  roughly  corresponds 
to  water  particles  in  air.  As  S  decreases,  the  solution  becomes  an  exponentially  damped 
oscillation  (Fig.  lb),  until  S  =  1.75  ((3  =  0.6667)  at  which  time  the  damping  disappears  (Fig. 
lc).  As  S  is  further  decreased,  the  solution  becomes  an  exponentially  increasing  oscillation 
(Fig.  Id). 

In  this  example  the  Stokes  drag  law  is  used  even  though  the  particle  relative  Reynolds  number 
greatly  exceeds  one.  With  Stokes  law  the  drag  is  essentially  zero  for  large  particle  Reynolds 
numbers.  Because  the  drag  is  essentialy  zero,  the  Basset  term  dominates  the  equation  of  motion. 
In  actuality  the  particle  drag  is  not  zero  in  the  region  with  high  particle  relative  Reynolds 
numbers  when  compressibility  and  inertia  effects  are  included  and  when  the  effect  of  the  Basset 
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term  is  reduced.  Forney  et  al.  (Ref.  14)  demonstrated  that  the  Basset  term  has  little  effect 
near  the  step  but  causes  the  particle  velocity  to  relax  more  slowly  away  from  the  step,  greatly 
increasing  the  particle  relaxation  distance.  This  example  indicates  that  the  Basset,  term  cannot 
be  arbitrarily  ignored  in  flows  with  large  velocity  gradients  even  when  fairly  dense  particles 
are  considered. 

A  second  example  of  the  importance  of  various  terms  in  Tchen’s  equation  of  motion  can 
be  seen  in  the  solution  given  by  Hjelmfelt  and  Mockros  (Ref.  2)  for  a  particle  responding 
to  a  uniformly  sinusoidally  oscillating  fluid  flow  field.  Their  results  for  amplitude  ratio  and 
phase  angle  of  the  particle  velocity  in  terms  of  the  fluid  velocity,  for  various  approximations 
of  the  equation  of  motion,  are 

—  =  J  0  +  fi)2  +  f22  (13) 

u  V 


and 


<f>  =  tan-1 


*2 


I  1 


Mi 


(14) 


where 


1.  General  solution 


1  + 


9  Nf 


1 


V2  (S  +  1/2)  I  S  +  1/2 


- — -  [zNf2  +  -^-1  2  +  fl  + 

(S  +  i/2)2  [  y/2  j 


9  Nf 


VT(3  +  1/2) 


(15) 


9  (L^)_  2  Nf2  +  Nf 


h  = 


(S  +  1/2)2 


V2 


81 


9  Nf 


-  |2Nf2  +  -^-1  +  1  +  - 

(S  +  1/2)2  I  V2"J  [  V2(S  +  1/2) 


(16) 


where 


S  =  gp/gf 


(17) 
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and 


Nf  = 


“V^ 


2.  Type  I  approximation,  neglecting  the  history  term, 


(18) 


f,  = 


(1  ~  S) 
(S  +  1/2) 


18Nf2 


(S  +  1/2) 


+  1 


18  Nf2 


(1  ~  S) 
(S  +  1/2)2 


ISNf2 


(S  +  1/2) 


-f  1 


3,  Type  II  approximation,  neglecting  the  history  and  apparent  mass  terms, 

f  =  (1  -  S)/S 


18Nf2 


+  1 


(S  +  1/2) _ 

ISNf2  (1  -  S)/S2 
ISNf2 


(S  +  1/2) 


+  1 


(19) 


(20) 


(21) 


(22) 


4.  Type  III  approximation,  neglecting  the  pressure,  history,  and  apparent  mass  terms, 

-1 


fi  = 


h  = 


18Nf2 


(23) 


S  J 

-ISN^/S 


+  1 


ISNf2 


(24) 


+  1 


Figure  2a  shows  the  amplitude  ratio  as  a  function  of  nondimensional  frequency  for  a  density 
ratio,  S,  of  1000,  which  approximately  represents  water  particles  in  air.  This  indicates  that 
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for  high  frequencies  or  for  large  particles  (small  Nf),  the  history  term  becomes  significant, 
whereas  the  apparent  mass  and  pressure  terms  are  insignificant.  The  phase  angle,  Fig.  2b, 
is  slightly  more  sensitive  to  the  additonal  terms.  This  example  corresponds  to  the  case  where 
the  relative  velocity  is  nearly  zero,  and  the  relative  acceleration  remains  finite.  As  the  density 
ratio  decreases,  as  indicated  in  Figs.  3a  and  b  for  S  =  2.65,  the  higher-order  terms  become 
more  significant. 

These  examples  indicate  that  the  Basset  term  cannot  be  ignored  in  a  general  sense  because 
of  order-or-magnitude  analysis  involving  only  the  coefficient  of  the  term,  as  is  usually  done. 
The  term  must  be  evaluated  on  a  case-by-case  basis  to  ensure  that  the  physics  of  the  particle/fluid 
interaction  are  properly  modeled.  This  is  especially  true  for  the  small,  light  particles  of  interest 
to  velocimetry  since  the  inertia  term  can  be  quite  small,  and  the  Basset  term  is  still  finite. 


2.3  EXTENSION  BEYOND  STOKES  FIDW 


The  repercussions  of  the  Stokes  flow  assumption  that  underlie  the  derivation  of  Eqs.  (1) 
and  (4)  should  be  noted.  The  assumption  is  that  the  convective  terms  of  the  Navier  Stokes 
equation  are  negligible.  This  assumption  limits  the  equation  to  very  low  Reynolds  numbers 
and  thus  to  cases  of  very  small  particle  lag.  Another  consequence  of  the  Stokes  flow  assumption, 
as  noted  by  Ho  (Ref.  15),  is  that  the  possibility  of  any  transverse  (lift)  forces  on  the  particle 
is  eliminated.  Saffman  (Ref.  16)  derived  the  lift  force  on  a  spherical  particle  in  a  uniform 
shear  flow  as 


Fl  =  1.6125  d2 


<3u 

dy 


i/2 

(v  -  u) 


(25) 


Evaluating  the  order-of-magnitude  of  this  force  relative  to  the  left-hand  side  of  Eq.  (4), 


Lift  Force 


'v  0 


-s 

A  _£f 
,d  Qp 


The  lift  force  is  then  of  the  same  magnitude  as  the  Basset  history  integral  for  a  particle  in 
a  shear.  Unfortunately,  there  is  no  general  expression  for  the  lift  force  on  a  sphere  in  an  unsteady 
nonuniform  flow  field,  and  hence  this  force  is  usually  ignored.  Fortunately,  the  lift  force, 
like  drag,  is  proportional  to  the  relative  velocity  and  is  usually  negligible  when  compared  to 
the  drag  force. 

The  general  procedure  for  extending  this  equation  to  higher  Reynolds  numbers  has  been 
to  multiply  various  terms  by  empirically  determined  coefficients.  Odar  and  Hamilton  (Ref. 
17)  performed  an  elaborate  experimental  study  and  concluded  that  empirical  constants  were 
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required  on  the  drag,  Basset,  and  apparent  mass  terms.  Later  work  by  Karanfilian  and  Kotas 
(Ref.  18)  indicated  that,  at  least  for  heavy  particles  in  sinusoidally  oscillating  flows,  the  drag 
term  required  empirical  modification,  and  the  other  terms  could  be  neglected.  In  this  work 
the  empirical  constants  of  Odar  and  Hamilton  (Ref.  17)  will  be  employed  because  of  the  more 
general  validation  and  more  thorough  scrutiny  they  have  been  given.  The  drag  term  in  the 
equation  of  motion  then  becomes 

Drag  Term  =  — m-  Cd  Wi|w,|  (26) 

4d 


where 


Wi  =  Vi  -  Uj  - 


24 


V2Uj 


and 


|wj|  =  S  Wj2 


The  empirical  coefficient  for  the  Basset  term  is  taken  from  Ref.  14  as 


Ah  =  0.48  + 


0.52  An 
(1  +  An)3 


(27) 


(28) 


(29) 


where 


d  d|wj 

|wj2|  dt 


(30) 


A  number  of  empirical  drag  coefficients  are  available  to  extend  the  particle  equation  of 
motion  beyond  the  Stokes  regime.  The  simplest  extensions  of  Stokes  drag  law  attempt  to  model 
the  incompressible  sphere  drag  coefficient  curve  (Fig.  4)  at  higher  Reynolds  numbers.  The 
more  complex  laws  add  corrections  for  compressibility,  rarefaction,  heat-transfer,  and  shape 
effects.  Donald  O.  Barnett  of  the  University  of  Alabama  presents  a  drag  coefficient  based 
on  the  experiments  of  Bailey  and  Hiatt  (Ref.  19)  of  the  form 

Cd  =  Cdo  .E  Fi  (31) 
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where 

Cd0  =  Stokes  drag  =  24/ReR 
Fi  =  Inertial  (Re)  effect 
F2  =  Compressibility  (M)  effect 
F3  =  Rarefaction  (Kn)  effect 
F4  =  Particle  shape  effects. 

Barnett  expresses  the  inertial  influence  coefficient  as 

Fj  (Rcr)  =  1  +  0.197  ReR  063  +  0.00026  ReR  (32) 

The  effect  of  increasing  Reynolds  number  is  to  increase  the  drag  above  that  which  is  predicted 
for  Stokes  flow,  as  seen  in  Fig.  5.  Compressibility  effects  are  modeled  by 


F2  (Mr)  =1+2/3  (Mr  -  0.15)2  (33) 

for  M  >  0.15,  As  shown  in  Fig.  6,  the  effect  of  compressibility  is  to  increase  the  drag  above 
Stokes  drag.  The  rarefaction  influence  coefficient  is  given  as 

F3  (Kn)  =  (1  +  Kn  (1  +  A  e-B/Kn))_i  (34) 

where  A  =  1.34285  and  B  =  0.277  for  Kn  <  0.85,  and  A  =  1.2243  and  B  =  0.18  for  Kn 
>  0.85.  As  shown  in  Fig.  7,  the  effect  of  rarefaction  is  to  decrease  the  Stokes  drag.  The  particle 
shape  is  accounted  for  by 


F4  (0)  =  1  +  4.065  (1  -  \p)Res  °-238 


where  \j/  is  the  shape  factor  defined  by 


t  = 


(35) 


(36) 


and 
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Ap  and  vp  are  the  surface  area  and  volume  of  the  particle.  Res  is  the  particle  Reynolds  number 
based  on  the  equivalent  spherical  diameter  from  Eq.  (37).  As  the  particle  becomes  less  spherical 
(or  as  \p  decreases  from  1),  the  drag  increases  as  shown  in  Fig.  8.  Heat-transfer  effects  may 
be  added  to  Barnett’s  drag  law  by  calculating  the  effective  viscosity  for  calculating  the  Reynolds 
number  as  recommended  by  Golovin  and  Fominysk  (Ref.  20).  This  is  done  by  calculating 
an  effective  temperature  that  is  the  geometric  mean  of  the  fluid  and  particle  temperatures. 
The  effective  viscosity  is  calculated  based  on  this  effective  temperature.  As  shown  in  Fig.  9, 
the  effect  of  increasing  the  temperature  of  the  particle  relative  to  the  temperature  of  the  fluid 
is  to  increase  sphere  drag. 

Barnett’s  drag  law  is  limited  to  particle  relative  Mach  number  (MR)  of  less  than  0.8, 
because  the  transonic  drag  rise  and  subsequent  supersonic  drag  of  a  sphere  is  not  modeled. 
Barnett’s  drag  law  is  useful  for  evaluating  the  effects  of  various  parameters  on  particle  motion 
but  must  be  used  with  caution  in  regions  of  large  particle  lag,  such  as  the  region  downstream 
of  a  strong  shock.  In  such  regions  the  drag  is  better  represented  by  the  drag  law  of  Crowe 
(Ref.  21). 

Cd  =  (Cdl  -  2)  exp  [-3.07  71/2  MR/ReR)G(ReR)]  (38) 

+  [h(MR)/71/2  Mr)]  exp  [  —  ReR/ (2MR)]  +  2 


where 


Cd,  -  (24/ReR)  (1  +  0.158  ReR2/3) 


(39) 


and 


log10  G(ReR)  =  1.25  [1  +  tanh  (0.77  log10  ReR  -  1.92)]  (40) 


and 


h(MR)  =  (2.3  +  1.7  (Tp/Tf)1/2)  -  2.3  tanh  (1.17  log)0  MR)  (41) 

Crowe’s  drag  law  includes  inertia,  compressibility,  rarefaction,  and  heat-transfer  effects. 
Although  Crowe’s  drag  law  is  valid  over  a  wider  range  of  particle  conditions  than  Barnett’s, 
it  makes  simple  parametric  studies  of  the  various  correction  terms  more  difficult. 
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2.4  NUMERICAL  SOLUTION  OF  THE  EQUATION  OF  MOTION 


Equation  (4)  can  be  put  into  the  following  form  to  facilitate  numerical  solution: 

3Q  mf 


dvj 

dt 


4d  (mp  +  1/2  mf) 


(mp  -  mf) 

Wi  W  4-  — — ~ - - —  D; 

(mp  +  l/2mf) 


mf 


Dui 


+ 


mf 


d 


(mp  +  l/2mf)  Dt  2(mp  +  l/2mf  dt 


Ui  [  Y  (t),t] 


d2  , 
+  —  V2 


40 


ui|Y(t)| 


Ah 


‘i'ird2n 


2(mp  +  I/2mf)  b 


}dT 

|  y/lTV(t  -  t)  J 


(42) 


This  can  be  interpreted  as  a  nonlinear  first-order  differential  equation  of  the  form 

dv; 


dt 


=  fi 


(43) 


A  few  comments  on  the  stability  of  Eq.  (42)  are  in  order.  If  we  assume  that  the  drag  is 
given  by  Stokes  law  for  simplicity,  and  we  ignore  the  lower  terms,  f;  can  be  written 


fi  = 


18  g 


d2  (6p  +  1/2  Gf) 


Wi 


The  Jacobian  is  then  given  by 


J  =  N, 


1  0  0 
0  1  0 
0  0  1 


where 


Ns  = 


18g 

d2  (ep  +  i/2gf) 


(4*4) 


(45) 


(46) 


18 


AEDC-TR-86-41 


The  system  of  equations  is  quite  stable  and  easy  to  solve  for  large  particles.  As  the  particle 
diameter  becomes  small,  N$  becomes  quite  large  as  shown  in  Fig.  10,  and  the  system  of 
equations  becomes  stiff.  Physically,  as  the  particle  diameter  approaches  zero,  the  relative  velocity 
also  approaches  zero.  Numerically  this  is  not  always  the  case,  and  hence  small  time  steps  are 
required  to  maintain  the  stability  of  the  system.  This  problem  can  be  circumvented  by  arbitrarily 
setting  the  particle  velocity  to  the  fluid  velocity  for  particle  sizes  below  a  given  size  or  for 
particle  relative  Reynolds  numbers  below  a  specified  value.  Since  interest  in  this  effort  is  in 
the  response  of  very  small  particles  to  turbulence,  a  limiting  value  of  particle  Reynolds  number 
of  10-12  was  chosen  to  be  compatible  with  the  convergence  criteria  used  in  numerically  solving 
the  particle  equation  of  motion. 


The  Basset  term  contains  an  apparent  singularity  at  the  upper  limit  of  integration 
(r  =  t).  This  must  be  considered  in  the  numerical  solution  of  the  particle  of  motion.  Setting 


f(r)  = 


(47) 


the  integral  is 

t  f(r)  dr 

O  yjt  —  T 


(48) 


Integrating  by  parts,  this  becomes 

2  VTf  (0)  +  2  \  VT^T-^-dr 


(49) 


When  the  integral  in  Eq.  (49)  is  evaluated  numerically,  the  contribution  to  the  sum  by  the 
point  r  is  zero.  Thus  this  point  is  not  a  singularity.  Physically,  this  is  consistent  with  the  proper 
interpretation  by  the  Basset  history  term  since  the  contribution  of  the  present  time  to  the 
history  of  the  particle  should  be  neglibie.  This  simplifies  the  numerical  solution  of  the  particle 
equation  of  motion  since  numerical  calculation  of  the  integral  term  does  not  require  a 
knowledge  of  the  velocities  at  the  present  time  step. 

A  computer  code  was  written  to  solve  Eq.  (42).  The  code  uses  a  fourth-order  forward 
marching  Adam’s  predictor-corrector  technique  to  integrate  for  the  particle  velocity  at  a  given 
time  step.  The  corrector  term  was  iterated  until  the  desired  convergence  was  obtained.  The 
convergence  criteria  were 

---  —  V'  *  <  1  x  10~6  (50) 

u1 
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where  i  was  the  iteration  index.  Each  component  of  the  velocity  was  required  to  converge 
independently.  These  criteria,  though  somewhat  restrictive,  were  found  to  guarantee  converged 
solutions  for  any  case  to  which  they  were  applied.  The  same  numerical  technique  was  used 
to  integrate  the  particle  velocity  to  obtain  particle  position.  The  integral  term  was  evaluated 
using  a  trapezoidal  rule  integration,  and  the  required  fluid  derivatives  were  evaluated  using 
simple  backward  differencing.  The  Adam’s  method  does  suffer  occasional  numerical  instability 
because  of  the  stiffness  of  the  equation  of  motion.  To  overcome  this  problem,  a  fourth-order 
Runge-Kutta  solver  was  used  to  start  the  solution  and  to  restart  the  solution  when  the  Adam’s 
method  diverged.  The  Adam’s  method  is  preferred  to  the  Runge-Kutta,  because  it  is  much 
faster  and  allows  simple  testing  for  convergence  to  a  desired  accuracy. 

3.0  PARTICLE  SIZE  DISTRIBUTION  EFFECTS 

A  commonly  occurring  and  often  ignored  phenomenon  in  laser  velocimetry  is  the 
broadening  of  measured  velocity  distributions  attributable  to  a  particle  size  distribution  in 
the  presence  of  a  velocity  gradient.  This  broadening  is  often  interpreted  as  an  increase  in 
turbulence  intensity  when  it  actually  has  nothing  to  do  with  the  turbulence  level  of  the  fluid. 
Three  simple  examples  frequently  found  in  high-speed  LV  applications  are  given  here. 

Before  proceeding  to  a  discussion  of  particle  size  distribution  effects  on  turbulence 
measurements  with  the  LV,  a  brief  outline  of  the  interpretation  of  velocimeter  data  will  be 
given.  The  discussion  will  be  limited  to  two-dimensional  data  for  simplicity,  although  extension 
to  three  dimensions  is  straightforward. 

A  sample  velocity  correlogram  is  given  in  Fig.  11.  For  typical  velocity  distributions,  the 
probability  distribution  takes  the  form  of  an  ellipse  when  plotted  in  velocity  space.  The 
turbulence  intensity  in  the  x  and  y  directions  is  given  by  the  variance  of  the  velocity  distribution 
in  each  direction  and  may  be  interpreted  as  the  square  of  half  the  projected  ellipse  length 

on  each  axis.  The  Reynolds  shear  is  calculated  from  the  cross  correlation  of  the  simultaneously 

measured  ux  and  uy  velocity  components  and  may  be  interpreted  as  the  angular  rotation  of 
the  correlogram  relative  to  the  coordinate  system  of  interest. 

The  determination  of  the  principle  turbulence  intensities,  or  those  where  the  Reynolds 
shear  is  zero,  is  given  by  the  following  relations  in  two  dimensions: 

uio  =  uj2  cos 26  +  U1U2  sin  26  +  U22  sin20  (51) 

U20  =  rq2  sin20  -  U1U2  sin  26  +  u22  cos20  (52) 
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where 


0  =  — -  tan-1 
2 


2  uju2 
ui2  +u2'2 


(53) 


and  u10'2  and  u20'2  are  the  principle  stresses  and  ui'2,  u2'2,  and  ui'u2'  are  the  stresses  in  the 
measurement  coordinate  system. 


3.1  NORMAL  SHOCKS 


A  very  simple  example  of  how  particle  dynamics  can  affect  turbulence  measurements  with 
the  LV  is  a  normal  shock.  The  momentum  of  a  particle  cannot  be  dissipated  in  the  few  mean 
free  paths  comprising  the  shock  thickness.  Theoretical  particle  response  to  a  bow  shock  in 
front  of  a  cylinder  at  a  Mach  number  of  four  is  shown  in  Fig.  12  for  olive  oil  particles  of 
various  sizes.  The  different  particle  sizes  respond  to  the  shock  at  different  rates.  The  velocity 
distribution  measured  using  a  polydispersed  particle  size  distribution  would  initially  be  tight 
(the  variance  would  be  small)  in  the  free  stream  ahead  of  the  shock.  The  distribution  rapidly 
spreads  as  the  particles  relax  at  different  rates  behind  the  shock.  Finally,  as  the  particles  begin 
to  relax  to  the  velocity  behind  the  shock,  the  distribution  forms  a  peak  at  the  low  velocity 
end  with  a  tail  to  the  higher  velocities.  As  the  larger  particles  relax  to  the  fluid  velocity,  this 
tail  begins  to  disappear,  and  the  low  velocity  peak  grows. 


Experimental  results  for  this  case  are  shown  in  Fig.  13.  The  flow  was  seeded  with  olive 
oil  using  a  Collison  atomizer,  which  produces  a  broad  size  range  of  particles  as  shown  in 
Fig.  14.  The  velocity  distributions  correspond  to  what  one  would  expect  from  a  particle 
dynamics  point  of  view.  There  is  a  strong  peak  upstream  of  the  shock  in  the  free  stream, 
a  very  broad  velocity  distribution  behind  the  shock,  and  a  velocity  peak  that  begins  to  develop 
at  the  lower  velocities  as  the  cylinder  is  approached.  The  actual  fluid  turbulence  values  are 
masked  by  the  spread  in  the  velocity  distribution  attributable  to  the  particle  size  distribution 
present  in  the  experiment. 

3.2  OBLIQUE  SHOCKS 


The  previous  example  can  be  extended  to  two-dimensional  flows  by  examining  the  simple 
case  of  an  oblique  shock.  Calculations  for  a  40-deg  oblique  shock  with  an  upstream  Mach 
number  of  2.2  are  shown  in  Figs.  15  and  16.  Xs  is  the  distance  from  the  shock  measured 
perpendicular  to  the  shock.  The  calculations  show  a  behavior  for  each  velocity  component 
similar  to  that  seen  for  the  case  of  the  normal  shock.  The  velocity  distribution  is  tight  and 
near  the  upstream  velocity  near  the  shock,  then  spreads  dramatically,  and  finally  develops 
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a  tight  peak  at  the  downstream  velocity  with  an  associated  tail  of  lagging  particles. 

Figures  16a  to  16e,  presented  in  velocity  space,  indicate  that  a  particle  size  distribution  will 
also  affect  the  Reynolds  shear  measurement.  As  mentioned  in  the  introduction  to  this  section, 
the  orientation  of  the  turbulence  ellipse  in  velocity  space  is  an  indication  of  the  Reynolds 
shear.  These  figures  indicate  that  the  ellipse  would  be  rotated  as  though  a  strong  negative 
Reynolds  shear  were  present.  This  particle-induced  velocity  correlation  could  well  hide  the 
actual  turbulence-induced  velocity  correlation. 

3.3  LAMINAR  SUPERSONIC  BOUNDARY  LAYER 

The  previous  two  cases  indicate  that  particle  dynamics  can  greatly  affect  turbulence 
measurements  with  the  LV  in  the  presence  of  strong  velocity  gradients,  such  as  the  gradients 
associated  with  shocks.  The  present  example  illustrates  that  this  effect  can  also  be  seen  in 
flows  with  milder  gradients,  such  as  in  the  laminar  boundary  layer  on  a  7-deg  half-angle  cone 
at  a  Mach  number  of  four. 

Figure  17  shows  calculated  particle  response  and  actual  LV  data  taken  in  a  laminar  boundary 
layer  on  a  cone  with  a  similar  seed  particle  size  distribution  as  shown  in  Fig.  14.  The  calculated 
results  in  Fig.  17  indicate  that  the  apparent  increase  in  turbulence  in  the  center  of  the  boundary 
layer,  noted  by  the  broadening  of  the  velocity  distribution,  is  actually  attributable  to  the  broad 
particle  size  distribution  present  and  is  not  a  feature  of  the  flow. 

3.4  MEAN  FLOW  PARTICLE  SIZE  REQUIREMENTS 

These  three  examples  point  out  the  care  that  must  be  taken  in  acquiring  and  interpreting 
LV  data.  Particle  dynamics  effects  must  be  considered  and  seed  particle  size  distributions  must 
be  minimized  if  meaningful  results  are  to  be  acquired  with  the  LV.  It  also  points  out  that 
the  LV  may  not  be  the  instrument  of  choice  for  making  turbulence  measurements  in  shock 
interaction  regions,  an  application  of  velocimetry  many  seem  to  make.  In  these  regions  care 
must  be  taken  to  identify  the  actual  particle  size  distribution  the  LV  is  “seeing”  and  to  assure 
that  measurements  are  made  outside  of  the  relaxation  region  for  these  particles.  These  examples 
do  not  indicate  what  size  particle  is  required  to  obtain  meaningful  turbulence  measurements 
since,  to  this  point,  only  mean  flow  particle  lag  has  been  considered. 

A  simple  expression  for  estimating  mean  flow  particle  lag  can  be  derived  from  the  particle 
equation  of  motion,  assuming  that  only  the  drag  term  is  important.  This  yields  the  following 
equation: 
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ui  _  ,  ,  cdQ  1  dvj 

Vi  Cd  NsVj  dt 

where  the  derivative  can  be  expressed 

dvj  _  dvj 

dt  Vj  dxj 


(54) 


(55) 


and  Cd0  is  the  Stokes  drag.  This  equation  is  not  considered  accurate  enough  for  data 
correction  but  is  valuable  for  identifying  areas  in  the  flow  in  which  particle  lag  may  be  present. 
The  equation  demonstrates  the  importance  of  maximizing  the  Stokes  number  in  regions  of 
the  flow  where  large  velocity  gradients  are  present. 


4.0  TURBULENCE  MODEL 


In  order  to  further  study  the  effects  of  particle  size  on  their  response  to  turbulence,  a 
model  of  the  turbulent  flow  must  be  constructed.  There  are  basically  two  ways  to  visualize 
turbulence  in  a  flow  field,  either  from  an  Eulerian  or  a  Lagrangian  reference  frame.  If  we 
were  to  visualize  turbulence  as  a  finite  number  of  discrete  eddies,  each  with  its  own  random 
fluid  characteristics  convected  with  mean  velocity  u,  then  in  one  dimension  the  turbulence 
may  be  represented  as  in  Fig.  18a.  As  these  eddies  are  convected  past  a  probe,  the  velocity 
measured  would  be  as  shown  in  Fig.  18b.  This  could  easily  be  represented  as  a  Fourier  series 
and  could  be  simulated  using  white  noise  Monte  Carlo  techniques.  Thus,  in  an  Eulerian 
reference  frame,  the  turbulence  appears  to  be  the  superposition  of  random  frequency  sine 
waves.  If  the  probe  were  given  a  velocity  relative  to  the  turbulence,  the  frequency  domain 
of  the  Eulerian  turbulence  spectra  would  change,  producing  an  apparent  higher  frequency 
if  the  probe  relative  velocity  is  of  opposite  sign  to  the  mean  fluid  velocity  u.  Hence,  the  Eulerian 
representation  of  turbulence  is  a  function  of  the  rate  of  convection  of  eddies  past  the  probe 
and  the  individual  eddy  characteristics. 

Particles  interact  with  turbulence  from  a  Lagrangian  point  of  view.  Figure  19a  is  an  example 
of  how  the  eddies  appear  to  the  particle.  At  time  t  the  particles  and  the  eddy  are  both  centered 
at  point  c  with  a  relative  velocity  of  Wj.  As  time  increases,  the  particle  exits  the  eddy  and 
enters  a  second  eddy  with  a  relative  velocity  w2.  In  the  Lagrangian  reference  frame,  the 
turbulence  appears  more  as  a  randomly  spaced  set  of  velocity  discontinuities  as  the  particle 
moves  from  one  eddy  to  the  next,  as  shown  in  Fig.  19b.  The  number  of  these  discontinuities 
encountered  by  the  particle  is  a  function  of  the  time  and  length  scales  of  the  turbulence  as 
well  as  the  size  and  density  of  the  particle.  If  a  particle  is  small  and  light  enough,  it  may 
remain  in  the  same  eddy  throughout  the  eddy’s  lifetime  and  see  very  few  discontinuities  as 
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it  passes  through  the  flow  field.  If  a  particle  is  large  and  heavy,  it  will  “bust  through”  the 
eddies  and  not  have  time  to  respond  to  any  of  the  velocity  discontinuities  it  encounters. 

For  the  large,  heavy  particles  that  fail  to  respond  to  the  turbulence,  modeling  of  the  time 
scale  and  the  internal  structure  of  the  turbulent  eddy  are  unimportant,  because  the  particle 
will  generally  encounter  a  new  eddy  at  each  time  step.  Several  simple  turbulence-particle 
interaction  models  take  advantage  of  this  fact  since  it  simplifies  the  calculation  of  particle 
response.  As  the  particles  become  smaller  and  lighter,  the  internal  structure  of  the  eddy  and 
the  time  and  length  scales  become  more  important.  Since  information  on  these  length  scales 
and  on  the  internal  eddy  structure  is  difficult  to  come  by,  the  proper  simulation  of  the  particle 
response  for  the  small  particles  of  interest  to  velocimetry  becomes  more  difficult. 

The  turbulence  model  used  in  this  effort  assumes  that  the  turbulence  is  multinormal  with 
the  off-diagonal  correlation  functions  given  by  the  Reynolds  shear  values.  The  mean  fluid 
velocities,  fluid  properties,  turbulence  intensities,  and  Reynolds  shears  are  assumed  to  be  given 
functions  of  position  in  the  flow  field,  determined  from  experimental  data  or  from  calculations. 


Information  on  time  scales  and  eddy  internal  structure  is  often  more  difficult  to  come 
by.  For  wake  flows  or  free  shear  flows,  empirical  relationships  are  available  from  Shuen  (Ref. 
3)  and  Tennekes  and  Lumley  (Ref.  22).  The  time  scale  is  defined  by 


TEi  =  7 

O 


Ui(r)Ui(t  -  r) 


dr 


(56) 


and  the  length  scale  is  defined  by 


oo  (ur.)  x  =  A  (ur.)  x  =  B 
Lh  =  [  - =— 1 - dr 


U: 


(57) 


The  Taylor  microscale  is  defined  as 


U 

(3ui/3xj)2 


These  can  be  approximated  by 


Xi  =  86.06 

Vu 


(59) 
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L„  =  0.077  A  -  Xl2  — 


L22  =  L33  =  — -  Ln 
2 

t  _  ^  Ln 

El  V3  vu 


(60) 

(61) 

(62) 


where 


3  _ 

U  =  E  u'2  (63) 

i=l  i 

and  A  is  a  constant  of  order  one.  These  scales  can  be  approximated  from  eddy  viscosity  or 
mixing  length  information  using  the  following  relation  from  Tennekes  and  Lumley  (Ref.  22): 


(64) 


where  Ln  is  assumed  to  be  the  direction  aligned  with  the  mean  flow  vector,  and  l  is  Prandtl’s 
mixing  length.  With  scales  defined  as  above,  the  turbulent  eddies  are  modeled  as  follows. 
The  random  duration  of  the  eddy  is  modeled  by  assuming  that  the  autocorrelation  is  given 
by  the  Poisson  from  e~t/TEi.  The  spatial  distribution  of  the  turbulence  must  also  be  modeled 
since  the  solid  particle  will  not  exactly  follow  any  single  fluid  particle.  The  spatial  distribution 
is  assumed  to  be  given  by  Frenkiel’s  analytical  form  (Ref.  23), 


UjcUjx 


exp 


-r 


(n  +  l)Ljj 


cos 


nr 


(n2+l)Lij 


(65) 


where  r  is  the  distance  between  the  fluid  particle  at  point  c  and  the  solid  particle  at  point 
x.  Currently  n  is  set  to  zero.  The  velocity  components  at  points  c  and  x  are  assumed  to  be 
governed  by  a  multinormal  distribution  with  a  covariance  matrix  of 


Ucl'Ud'  Ud'Urf'  ud'uxl'  ucl'ux2' 

Uc2'Ucl'  UC2 '  UC2 '  Uc2 '  Uxl '  Uc2 '  Ux2 ' 
Uxl'Ud'  Uxl'uc2'  Uxl'uxl'  UX,'UX2' 
Ux2'ucl'  ux2 ' uc2 '  ux2'uxi'  ux2'ux2' 


(66) 
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The  turbulent  intensities  and  Reynolds  shear  at  point  x  are  assumed  to  be  the  same  as  those 
at  point  c  because  of  the  small  size  of  the  eddy  relative  to  the  scale  of  the  mean  flow.  The 
cross  terms  between  points  x  and  c  are  assumed  to  be  given  by  Eq.  (65).  The  random  fluid 
velocities  at  points  x  and  c  are  generated  by  the  multinormal  form  outlined  by  Rubinstein 
(Ref.  24), 


where 


Ulc 

U2c' 

Ulx' 

u2x' 


—  CijRzi 


Qj  = 


- j-i 

ui  uj  Mk'-'jk 


j-1 

E 

k  =  l 


(67) 


(68) 


Rzi  are  randomly  generated  one-dimensional  normal  variates.  The  correlation  given  by  Eq. 
(64)  begins  to  lose  meaning  at  large  values  of  r,  and  the  random  velocities  at  point  x  become 
independent  of  the  velocities  at  point  c.  To  avoid  losing  this  model  of  the  internal  structure 
of  the  turbulence,  the  code  changes  from  following  the  fluid  particle  at  point  c  to  following 
the  fluid  particle  at  point  x  if  r  is  greater  than  the  Taylor  microscale  Ai. 


The  variates  Rzi  and  Rz2  are  generated  when  an  eddy  is  created,  and  they  remain  constant 
for  the  lifetime  of  the  eddy.  The  variates  R23  and  RZ4  are  generated  at  each  time  step.  If  the 
eddy’s  random  lifetime  is  exceeded,  then  new  random  variates  Rzi  and  Rz2  are  generated,  and 
the  model  begins  to  follow  the  fluid  particle  located  at  point  x.  If  the  distance  r  between  points 
x  and  c  exceeds  the  dissipation  length  scale  X,  then  the  model  changes  from  following  the 
fluid  particle  at  point  c  to  the  fluid  particle  at  point  x.  This  is  accomplished  by  calculating 
new  random  variates  Rzi  and  Rz2  from  the  fluid  velocity  at  point  x  from 


Ulx 

Cn 


(69) 


Rz2  = 


u2x~  Cl2  Rzi 

C22 


(70) 
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These  variates  remain  constant  until  the  length  or  time  scales  of  the  eddy  are  exceeded. 

Figure  20  graphically  describes  the  process  the  model  uses  in  following  a  particle  through 
a  turbulent  flow  field.  In  Fig.  20a,  the  particle  begins  at  the  center  of  a  fluid  eddy.  The  minimum 
radius  of  this  eddy  is  the  Taylor  microscale  Xj.  As  time  progresses,  the  particle  remains  in 
this  eddy  until  the  eddy’s  random  life,  governed  by  a  Poisson  distribution,  is  exceeded.  The 
process  is  modeled  by  drawing  a  random  variate  from  [0,1]  at  each  time  step  and  comparing 
it  to  At/TEi.  If  the  variate  is  less  than  At/TE1,  new  values  of  Uj'  are  generated;  if  not,  Uj' 
remains  the  same  during  time  interval  At. 

In  Fig.  20b,  the  particle  again  begins  at  the  center  of  an  eddy  of  minimum  radius  X( .  As 
time  progresses,  the  particle  exits  this  eddy.  At  this  point,  the  code  begins  to  follow  the  eddy 
centered  at  point  x.  This  proceeds  until  the  particle  exits  the  new  eddy,  or  the  random  duration 
of  the  eddy  is  exceeded. 

The  turbulence  model  used  in  this  effort  attempts  to  model  the  internal  structure  of  the 
turbulence.  For  the  small,  light  particles  required  for  accurate  velocimetry  measurements,  it 
is  this  structure  that  the  particles  must  follow.  The  structure  is  somewhat  simplified  by 
assumptions  concerning  the  probability  distributions  governing  the  process.  Although  it  is 
fully  realized  that  all  turbulence  velocity  distributions  are  not  truly  normal,  and  length  and 
time  correlations  are  not  strictly  Poisson,  these  forms  should  be  adequate  to  describe  the  particle 
response  within  the  framework  of  a  simple  turbulence  model. 

5.0  PARTICLE  RESPONSE  TO  TURBULENT  FLOW  FIELDS 

The  particle  equation-of-motion  code  developed  in  Section  2.0  was  coupled  with  the  Monte 
Carlo  turbulence  model  developed  in  Section  4.0.  Individual  particles  were  followed  through 
the  randomly  varying  flow  field.  One  thousand  samples  were  used  to  calculate  the  particle 
and  fluid  statistics  at  various  stations  in  the  flow.  The  number  1000  was  chosen  as  a  compromise 
to  allow  a  reasonable  confidence  band  without  consuming  excessive  computational  resources. 

5.1  GRID-GENERATED  TURBULENCE 

The  code  was  first  applied  to  the  grid-generated  turbulence  experiment  of  Snyder  and 
Lumley  (Ref.  25).  This  provided  a  flow  field  with  near-homogeneous  and  near-isotropic 
turbulence  with  no  mean  flow  velocity  gradients.  In  this  experiment  particle  velocity 
autocorrelation  functions  and  particle  dispersion  were  measured  in  a  vertical  6.55-m/sec  flow 
behind  a  2.54-cm-sq  mesh.  Fluid  turbulence  intensities  were  measured  with  a  hot  wire  and 
are  shown  in  Fig.  21.  These  curves  are  shown  to  be  adequately  represented  by 
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42.4  (x/m  -  16) 


(71) 


39.4  (x/m  -  12) 


(72) 


where  m  is  the  mesh  size.  The  Lagrangian  integral  length  scale  and  Taylor  microscale  were 
tabulated  in  Ref.  25.  These  tabulated  values  were  fitted  by 

Ln  =  2.19327  +  0.0128493  x/m  +  8.582  x  10~6  2  (73) 

\t  =  0.00436  x/m  +  0.2507  (74) 


Equations  (61)  and  (62)  were  used  to  obtain  L22  and  TEi.  The  length  and  time  scales  are  shown 
in  Figs.  22  and  23,  respectively. 


Experimental  results  are  given  in  Ref.  25  for  the  four  different  particles  given  in  Table 
1.  The  Stokes  number,  defined  as, 

Ns  =  - ^ -  (75) 

d^  (gp  4-  l/2gf) 

is  given  as  a  parametric  indication  of  particle  response.  As  the  Stokes  number  decreases,  particle 
response  to  turbulence  will  also  decrease. 

Calculations  were  made  for  the  four  different  particles  in  Table  1.  The  particle  calculations 
were  begun  at  the  experimental  particle  injection  point  of  x/m  =  20.  The  particle  velocity 
was  assumed  equal  to  the  randomly  fluctuating  fluid  velocity  at  the  injection  point.  The 
experimental  measurements  began  at  x/m  =  41.  Because  of  the  large  separation  between  the 
injection  point  and  the  measurement  stations,  there  was  no  influence  of  the  assumed  particle 
initial  conditions  on  the  calculated  particle  response  at  the  measured  stations. 

Comparison  of  the  predicted  and  measured  particle  dispersions  is  shown  in  Fig.  24.  The 
agreement  in  all  cases  is  considered  good.  Even  better  agreement  can  be  achieved  for  each 
case  by  slightly  varying  the  Lagrangian  length  and  time  scales  for  each  different  kind  of  particle. 
This  was  not  done  because  such  “fine  tuning”  cannot  be  done  when  the  code  is  used  to  predict 
the  response  of  the  smaller  particles  of  interest  to  velocimetry.  Hence,  this  example  should 
be  representative  of  the  results  to  be  expected  from  this  code  for  various  size  and  density 
particles. 
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Predicted  and  experimental  turbulent  velocity  decays  are  shown  in  Fig.  25.  Again  there 
is  excellent  agreement.  The  roll  off  of  the  copper  particle  decay  predicted  by  the  code  is  not 
shown  by  the  experimental  data,  but  Snyder  (Ref.  25)  indicates  he  has  little  faith  in  the 
measurements  at  this  point  because  of  his  measurement  technique  and  the  very  low  values 
of  particle  velocity  fluctuations  there.  For  the  smaller,  lighter  particles  of  interest  in  this  study 
the  code  does  exceptionally  well. 

Following  the  successful  validation  of  the  code,  calculations  were  made  for  various  particle 
sizes  and  particle  densities  to  investigate  what  particles  were  required  to  properly  characterize 
the  turbulence.  Results  of  these  calculations  are  shown  in  Fig.  26  for  two  different  axial  stations. 
The  figure  indicates  that  the  particle  response  is  essentially  the  same  for  these  two  locations, 
which  is  to  be  expected  because  of  the  small  change  in  turbulent  intensities  between  them 
(a  turbulent  intensity  of  1.9  percent  at  x/m  of  80  compared  with  a  turbulent  intensity  of  1.5 
percent  at  x/m  of  124).  The  turbulence  is  properly  reproduced  by  particles  with  Stokes  numbers 
greater  than  1.3  x  104,  which  corresponds  to  a  5-/im  water  particle.  The  mean  velocity  values 
are  still  quite  good  for  particles  with  Stokes  numbers  of  12,  which  corresponds  to  a  160-^m 
water  particle.  The  calculations  also  indicate  that  the  particles  respond  differently  to  different 
components  of  the  turbulence.  The  particles  respond  to  the  turbulence  component  in  the 
direction  of  the  mean  flow  vector  significantly  better  than  they  respond  to  the  component 
of  turbulence  normal  to  the  mean  flow  vector.  This  difference  in  response  seems  to  disappear 
as  the  particles  begin  to  follow  the  fluid  velocity  turbulent  excursions.  The  results  indicate 
that  the  fluid  turbulence  properties  will  be  underpredicted  when  particles  that  are  too  large 
to  respond  to  the  flow  are  used. 

Some  physical  insight  into  the  particle  response  to  turbulence  can  be  gained  from  the 
Lagrangian  turbulence  described  in  Section  4.0  and  shown  in  Fig.  19.  The  particle  response 
is  indicative  of  the  ratio  of  the  particle  relaxation  time  to  the  time  scale  of  the  turbulence. 
If  the  particle  relaxation  time  is  of  the  same  order  as  the  time  scale  of  the  turbulence,  the 
particle  never  has  an  opportunity  to  relax  to  the  fluid  velocity  of  the  eddy  in  which  it  is  located. 
As  the  particle  relaxation  time  decreases,  the  particle  more  faithfully  reproduces  the  fluid 
motion.  The  Stokes  number  is  an  indicator  of  the  time  response  of  a  particle  and  serves  as 
a  good  parameter  for  evaluating  particle  response.  Figure  27  indicates  that  the  product  of 
the  Stokes  number  and  the  integral  time  scale  must  be  of  the  order  of  1000  for  accurate 
turbulence  measurements  in  this  flow  field. 
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5.2  SUBSONIC  AXISYMMETRIC  JET 

The  second  validation  case  for  the  code  was  the  axisymmetric  jet  experiment  of  Yuu  et 
al.  (Ref.  26).  In  this  experiment,  particle  diffusivity  was  measured  in  a  simple  self-preserving 
subsonic  jet.  This  flow  offered  the  opportunity  to  evaluate  the  effect  of  particle  dynamics 
on  Reynolds  shear  measurements  and  also  included  the  effects  of  a  mean  flow  velocity  gradient. 

Yuu  presents  curve  fits  of  the  turbulence  intensities  and  mean  flow  velocities  for  this  jet 
based  on  the  experiments  of  Liepman  and  Laufer  (Ref.  27),  Laurence  (Ref.  28),  Corrsin  and 
Uberoi  (Ref.  29),  and  Wygnanski  (Ref.  30).  The  basic  flow  structure  is  shown  in  Fig.  28.  The 
potential  core  of  the  jet  is  defined  by 


x/D  <6.8 


(76) 


r/D  <  0.5 - — 

13.6 


(77) 


where  D  is  the  nozzle  exit  diameter.  The  velocities  in  this  region  are  given  by 
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u0 
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u0 

“  0 


where  uQ  is  the  nozzle  exit  velocity.  In  the  mixing  region,  defined  by 
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the  velocities  are  given  by 
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—  =  -4.05  ^i2  -  11.7  r?,2  +  255  m4 


r2 

i  =  0.0217  exp  (-200  rj i2) 


(84) 

(85) 


2  =  0.0103  exp  (-217  rji2) 


(86) 


where  ijj  =  (f  —  0.5 )/x.  In  the  main  region  of  the  jet,  defined  by  x  >  6.8,  the  velocities 
are  given  by 


=  (6.8  -  630  r/22  +  2313  *?23)  x"1 


(87) 


u2 

u0 


_  =  (3.4  V2  ~  472  t,23  +  1851  r/24)  x 


v-i 


(88) 
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ui  =  — r-  exp  (-154  t]22) 


(89) 


u22  =  i^-exp  (-178  i?2) 


u, 


2  x2 


(90) 


where  =  r/x.  The  Reynolds  shear  in  the  mixing  region  was  taken  from  Chow  (Ref.  31)  as 


ui  u2  =  0.012  (1  -  204.5  vi2)  exp  [-(30  m  +  0.2)2/1.704] 


(91) 


Un 
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The  Reynolds  shear  in  the  main  jet  is  given  by  Hinze’s  expression  (Ref.  23)  as 


Ui  u2  _  23.12  r/2  r  _  r?22 

u02  V2  x2  [  0.01568 

The  integral  length  scale  is  approximated  by  (Refs.  26  and  30) 

(92) 

Lm  =  Lm/D  =  0.0385  x 

(93) 

Lh/d  =  Lm  (exp  j-100  f  +  0-45)/x]2 

) 

+  exp  |  -100  [(r  -  0.45)/x]2|)  1/2 

for  x  <  4, 

(94) 

Lu/d  =  Lm  exp  (-40  r/22) 

for  4  <  x  <  10, 

(95) 

Lu/d  =  Lm  exp  (-50  r) 22) 

For  x  >  10.  L22,  Xi,  and  Tm  are  given  by  (Ref.  30) 

(96) 

L22  =  Ln/2.54 

(97) 

The  above  fluid  velocity  and  turbulence  scales  are  graphically  represented  in  Figs.  29  through  34. 

Calculations  were  made  for  20-/zm-diam  fly-ash  particles  in  a  100-m/sec  jet  to  compare 
with  the  experimental  results  of  Yuu  (Ref.  26).  The  calculations  were  begun  with  the  particle 
in  the  potential  core  region  of  the  jet  moving  at  the  jet  fluid  velocity.  Predictions  of  the  mean 
particle  velocity  on  the  nozzle  centerline  are  compared  with  experimental  results  in  Fig.  35. 
As  shown,  the  agreement  is  quite  good.  Comparisons  of  predicted  and  measured  particle 
diffusivities,  defined  by 
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_  epi  =  y  ~  (tp  -  fp)2  (100) 

are  shown  in  Fig.  36.  Agreement  is  again  quite  good,  especially  considering  that  the  time 
derivative  of  a  statistically  determined  quantity  is  being  compared. 

Calculations  were  then  made  for  particles  in  the  size  range  of  interest  for  velocimetry. 
The  calculations  were  made  in  the  near  region  of  the  jet  (x/D  <  12.0)  to  minimize  the  time 
step  required  in  the  numerical  calculation  and  off  centerline  in  the  region  of  maximum  Reynolds 
shear.  Results  are  shown  in  Fig.  37.  The  general  trends  at  all  x  locations  are  the  same;  the 
mean  velocity  ratio  approaches  one  first  as  the  Stokes  number  increases,  followed  by  the  ratio 
of  the  longitudinal  turbulent  stresses,  the  ratio  of  the  lateral  turbulent  stresses,  and  finally 
the  ratio  of  the  Reynolds  shear.  These  calculations  are  for  5.0-,  2.5-,  1.0-,  0.25-,  and  0.1-/*m 
water  particles  in  regions  of  the  flow  where  the  turbulence  intensity  is  about  20  percent.  The 
particles  are  less  responsive  at  x/D  =  8.57,  because  the  magnitude  of  the  velocity  fluctuations 
is  approximately  50  percent  higher  there. 

The  effect  of  particle  dynamics  can  be  seen  in  the  velocity  correlograms  shown  in  Fig. 
38.  As  the  particle  size  increases,  the  ellipse  begins  to  shrink  and  to  rotate  so  that  the  principle 
axes  align  with  the  mean  flow  direction.  This  corresponds  to  an  underprediction  of  turbulence 
quantities,  which  worsens  as  the  particle  size  increases. 

The  calculations  indicate  that  a  water  particle  of  maximum  diameter  of  0.25  /tm  is  required 
to  make  reasonably  accurate  turbulence  measurements  in  this  flow  field  (reasonably  accurate 
being  defined  as  the  size  for  which  all  the  turbulence  quantity  particle  to  fluid  ratios  are  better 
than  98  percent),  whereas  a  2.5-^m  water  droplet  will  yield  acceptable  results  for  the  mean 
velocity.  Figure  39  indicates  that  the  product  of  the  Stokes  number  and  the  integral  time  scale 
must  be  greater  than  200  for  meaningful  turbulence  measurements  in  this  flow  field. 

5.3  MULTIMODAL  TURBULENCE 

The  turbulence  model  developed  here  is  based  on  an  analytical  representation  of  the  time- 
averaged  flow  field  that  is  consistent  with  a  single  normal  mode  representation  of  turbulence, 
but  it  is  not  always  consistent  with  actual  LV  measurements  of  turbulence.  This  can  be  seen 
in  the  two  stream-mixing  measurements  made  with  the  LV  reported  in  Ref.  1.  The  velocity 
correlogram  was  shown  to  be  bimodal  (dumbbell-shaped)  with  two  distinct  velocity  modes 
representing  the  high-speed  and  the  low-speed  streams.  The  two  modes  begin  to  coalesce  as 
the  mixing  region  more  fully  develops  and  eventually  become  one  mode  far  downstream.  In 
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the  early  region  of  the  mixing  zone  the  maximum  Reynolds  shear  is  very  large,  because  the 
stress  ellipse  has  a  very  long  major  axis  and  a  short  minor  axis.  This  is  true  even  if  the  individual 
components  of  the  dumbbell  are  isotropic.  If  the  particle  lag  effects  on  mean  velocity  are 
negligible  (as  was  shown  in  Ref.  1),  then  the  mean  velocity  of  each  component  of  the  dumbbell 
is  accurately  determined.  As  a  result  both  the  major  axis,  uio,  and  the  angle,  6,  of  the  stress 
ellipse  (Fig.  11)  are  accurately  determined.  Only  the  minor  axis,  U20,  of  the  stress  ellipse  is 
sensitive  to  particle  lag  and  would  be  underpredicted  if  the  particles  were  too  large  to  follow 
the  turbulence.  The  Reynolds  shear  in  the  mean  flow  direction  is 

-7-7  (U)O  -  U20) 

Uiu2  =  - — - sin  20  (101) 

and  will  be  overpredicted  when  particle  dynamics  effects  are  present. 

The  weighting  of  the  individual  components  of  the  dumbbell  was  shown  to  be  a  function 
of  the  respective  individually  controlled  seed  rate  present  in  each  stream  in  Ref.  1.  The  correct 
weighting,  or  intermittency,  of  the  velocity  modes  cannot  be  determined  from  LY  data  alone. 
To  properly  model  this  flow,  an  independent  measurement  of  the  fluid  intermittency  would 
be  required. 

6.0  APPROXIMATIONS  FOR  PARTICLE  RESPONSE  TO  TURBULENCE 

The  previous  results  indicate  that  quite  small  particles  may  be  required  to  obtain  accurate 
turbulence  measurements  with  the  LV  in  all  but  very  low-speed  flows.  It  would  be  useful  to 
know  what  size  particle  is  required  for  a  specific  flow  field  to  aid  in  pretest  planning  and 
posttest  analysis  of  LY  data.  The  numerical  method  described  in  Sections  4.0  and  5.0  is  not 
a  convenient  tool  for  general  velocimeter  applications  since  the  majority  of  the  required  flow- 
field  and  turbulence  information  is  not  known  a  priori.  Hence,  some  methods  of  approximating 
the  particle  response  to  turbulence  would  be  useful  in  practical  application.  Two  possible 
techniques  will  be  discussed  in  this  section. 

6.1  SPECTRAL  CUTOFF  FREQUENCY  APPROXIMATION 

The  work  of  Hjelmfelt  and  Mockros  (Ref.  2)  discussed  in  Section  2.0  offers  one  method 
of  approximation.  This  method  is  based  on  the  ability  of  a  particle  to  respond  to  the  highest 
frequency  present  in  a  turbulent  flow  field.  The  highest  frequency  in  the  Eulerian  time  spectrum 
can  be  approximated  by  (Ref.  22) 

1/2 

IT1  (102) 
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where 


r)  =  (v3/e)l/4 


The  dissipation,  e,  can  be  approximated  by  (Ref.  3) 

(0.09)3/2  U2 
40  v 


The  frequency  parameter  may  be  defined  as 


Nf0  =  4.315 


v 

d  VU 


(103) 


(104) 


(105) 


or,  in  terms  of  the  Stokes  number, 


Nf0  -  0.719 


v  Ns  (2S  +  1) 


U 


(106) 


The  particle  response  can  then  be  found  from  Eq.  (13)  for  various  size  and  weight  particles. 

6.2  SIMPLE  TURBULENCE  MODEL  APPROXIMATION 

Meek  and  Jones  (Ref.  32)  present  a  method  of  predicting  the  response  of  heavy  particles 
to  turbulence  based  on  a  simple  isotropic  turbulence.  This  simple  model  is  valid  when  the 
particles  are  much  more  dense  than  the  fluid  and  when  the  particles  are  small  enough  for 
Stokes  law  to  apply.  Their  equation  is 


where 


,'2 
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1  +  eu 


(107) 


en  = 


d2  (2S  +  1) 
180  Tei 


Ns  TEj 


(108) 


Substituting  into  Eq.  (62)  for  TE1,  en  becomes 


=  8.436  x  10-=  U  d2  <2S  +  S 

Av2 


3.037  x  10-3 


U 


ApN, 


(109) 
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The  particle  response  can  then  be  approximated  from  Eq.  (107)  for  various  size  and  weight 
particles.  This  equation  is  essentially  a  form  of  the  product  of  the  integral  time  scale  and 
the  Stokes  number,  which  forms  a  nondimensional  quantity  that  was  shown  to  correlate  the 
particle  response  results  well  in  Section  5.0. 

6.3  APPLICATIONS  OF  THE  APPROXIMATION  TECHNIQUES 

Both  techniques  were  applied  to  the  grid-generated  turbulence  and  the  axisymmetric  jet 
examples  discussed,  as  described  in  Section  5.0.  Results  for  the  grid  turbulence  are  shown 
in  Fig.  40;  jet  results  are  shown  in  Fig.  41.  These  results  indicate  that  the  method  of  Meek 
and  Jones  better  represents  the  predicted  trend  of  particle  response  with  decreasing  Stokes 
number.  The  method  of  Hjelmfelt  and  Mockros  predicts  a  much  greater  loss  of  response  with 
decreasing  Stokes  number  than  actually  exists.  Both  methods  yield  about  the  same  value  of 
Stokes  number  for  turbulent  intensity  ratios  of  0.99.  Both  methods  underpredict  the  correct 
minimum  value  for  the  grid  case;  whereas  they  both  predict  the  correct  minimum  Stokes  number 
for  the  jet.  For  the  cases  presented,  the  method  of  Meek  and  Jones  seems  to  be  adequate 
for  “rule  of  thumb”  approximations  of  the  expected  particle  response  to  turbulence  and  for 
estimating  LV  data  quality. 


7.0  SUMMARY 

A  numerical  model  for  evaluating  seed  particle  requirements  for  accurate  turbulence 
measurements  with  the  laser  velocimeter  has  been  developed.  The  particle  equation  of  motion 
was  carefully  analyzed  and  the  importance  of  including  the  Basset  history  integral  noted.  A 
Monte  Carlo  turbulence  model  was  developed,  which  included  an  elementary  description  of 
the  internal  structure  of  an  eddy.  When  these  models  were  applied  to  simple  flow  fields  to 
evaluate  particle  dynamics  effects  on  LV  turbulence  measurements,  two  important  conclusions 
were  reached. 

The  first  conclusion  stems  from  an  analysis  of  a  particle  size  distribution  present  in  a 
flow  field  with  sufficient  mean  velocity  gradients  to  cause  the  particles  to  lag  the  fluid  velocity. 
The  response  of  the  individual  particle  sizes  will  lead  to  a  broadening  of  the  measured  velocity 
distribution.  This  broadening  can  affect  both  the  turbulent  intensity  and  Reynolds  shear 
measurements.  These  particle  dynamics  effects  also  point  to  limitation  of  the  velocimeter  for 
obtaining  turbulence  information  in  the  region  immediately  behind  shocks  or  in  other  flows 
with  large  velocity  gradients. 

Secondly,  particles  that  are  capable  of  tracking  the  mean  flow  may  not  be  adequate  for 
measuring  turbulence  quantities.  Large  monodisperse  particles  will  cause  the  turbulence 
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quantities  to  be  underpredicted.  As  the  particle  diameter  or  density  decreases,  the  fluid  mean 
velocity  will  be  reproduced  first,  followed  by  the  longitudinal  turbulent  stress,  the  lateral 
turbulent  stress,  and  finally  the  Reynolds  shear.  It  is  anticipated  that  the  higher  order  the 
turbulence  correlation  parameter,  the  more  sensitive  the  measurement  will  become  to  particle 
dynamics. 

Calculations  were  made  for  two  simple  turbulent  flows,  grid-generated  turbulence  and 
a  subsonic  axisymmetric  jet.  The  Stokes  number  was  shown  to  be  a  reasonable  correlation 
parameter  for  particle  diameter  and  density  effects.  It  was  found  that  the  product  of  the  Stokes 
number  and  the  turbulence  integral  time  scale  (which  is  equivalent  to  the  ratio  of  the  turbulence 
integral  time  scale  to  the  particle  response  time)  must  be  of  order  1000  for  grid-generated 
turbulence  and  200  for  the  jet  turbulence  for  meaningful  turbulence  measurements  with  the 
LV.  This  leads  to  a  requirement  for  submicron  particles  for  most  applications  of  the  LV. 
Approximation  techniques  for  particle  response  to  turbulence  were  developed  for  “rule  of 
thumb”  type  applications. 

It  is  interesting  to  note  that  the  particle-induced  errors  associated  with  the  LV  have  a 
cancelling  effect.  The  mean  flow  particle  lag  leads  to  a  broadening  of  the  measured  velocity 
distribution  when  a  particle  size  distribution  is  present,  whereas  large  particles  will  under¬ 
respond  to  the  actual  fluid  turbulence.  In  complex  flow  fields,  both  of  these  effects  may  be 
occurring,  and  care  should  be  taken  to  assure  that  measurements  that  appear  reasonable  are 
not  the  result  of  a  fortuitous  interaction  between  these  competing  particle  dynamics  effects. 

Several  recommendations  for  future  work  in  this  area  can  be  made.  There  is  a  very  limited 
data  base  for  particle  response  to  turbulence.  Carefully  controlled  experiments  in  known,  well- 
documented  flow  fields  with  smaller  size  particles  would  aid  in  verifying  this  code  and  improving 
the  turbulence  model.  Numerical  improvements  in  the  particle  equation-of-motion  solver,  such 
as  the  use  of  a  stiff  system  equation  solver  and  faster  integration  techniques  in  the  evaluation 
of  the  Basset  history  integral,  would  reduce  the  cost  and  manual  intervention  required  for 
solutions  with  this  code. 

F.  L.  Heltsley’s  study  (Ref.  1)  included  bimodal  results  in  the  initial  region  of  a  jet/base 
and  a  free-stream/base  shear  layer.  The  bimodal  nature  of  the  velocity  distributions  is  quite 
pronounced  near  the  jet  exit  and  subsides  as  the  mixing  region  grows.  The  separate  velocity 
modes  form  a  “dumbbell”  when  plotted  in  velocity  space,  each  mode  containing  a  distinct 
mean  and  turbulence  level  consistent  with  fluid  streams  present  on  either  side  of  the  shear 
layer.  The  mean  velocity  and  turbulence  levels  associated  with  a  point  in  the  shear  layer  are 
generated  by  processing  the  velocity  data  for  both  modes,  giving  a  somewhat  distorted  view 
of  the  actual  fluid  behavior.  The  relative  size  of  the  modes  could  be  controlled  by  regulating 
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the  amount  of  seed  introduced  to  each  stream  individually.  This  creates  a  perplexing  situation 
since  the  mean  and  turbulence  measurements  with  the  LV  become  a  function  of  the  seed  rate 
from  the  two  seeders.  Further  work  is  required  to  assure  that  the  multimodal  measurements 
are  real  and  not  particle-induced  and  to  develop  data  sampling  techniques  for  obtaining  LV 
measurements  in  such  flows. 

Finally,  there  are  some  other  factors  not  discussed  in  this  work  that  may  also  effect 
turbulence  measurements  with  the  LV,  including  machine  (processor)  broadening,  preferential 
treatment  by  the  processor,  velocity  biasing,  and  spurious  noise  in  the  measurement  systems. 
A  detailed  study  of  these  effects,  coupled  with  this  work,  would  establish  the  accuracy  limits 
of  turbulence  measurements  with  the  LV  and  aid  in  the  interpretation  of  LV  data. 
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Velocity,  v/ 


a.  S  =  1000 

Figure  1.  Particle  response  to  a  step  change  in  fluid  velocity. 
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b.  S  =  5 

Figure  1.  Continued. 
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c.  S  =  1.75 
Figure  1.  Continued. 


AEDC-TR-B6-41 


Velocity,  v/a 


1.0  2.0  3.0  4.0  5.C 

O 

Time,  sec  x  10" 


d.  S  =  1.0 
Figure  1.  Concluded. 
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Figure  6.  Empirical  correction  for  subsonic  compressibility  effects. 


Rarefaction  Influence  Coefficient,  F,(Kn) 


Knudsen  Number,  Kn 

Figure  7.  Empirical  correction  for  rarefaction  effects. 
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Figure  13.  Measured  LV  velocity  probability  distributions  for 
a  cylinder  bow  shock  at  =  4.0  and  Ren  =  10s. 
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Figure  14.  Particle  size  distribution  for  olive  oil  in  a  Collison  Atomizer 
measured  with  an  ASAS-X®  Aerosol  sampling  unit. 


57 


AEDC-TR-86-41 


a.  x  Velocity 

Figure  15.  Particle  relaxation  behind  an  oblique  sbock  at 
Mg,  =  2.2  and  Be  =  4.4  X  10VM. 


b.  y  Velocity 
Figure  15.  Concluded. 
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a.  xM  =  0.015  mm 

Figure  16.  Particle  response  to  an  oblique  shock  for  Mi  =  2,2, 
0  =  40  deg,  6  =  14  deg,  and  Re  =  4.4  X  10VM. 
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Figure  16.  Concluded. 
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Eddies 


a.  ID  depiction  of  turbulence 


b.  TUrbuIent  velocity  as  seen  by  a  probe 
Figure  18.  Eulerian  representation  of  turbulence. 


Probe 
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a.  Turbulent  field  as  seen  by  a  small  particle 


t 

b.  Ttarbulent  fluid  velocity  as  seen  by  a  small  particle 
Figure  19.  Lagrangian  turbulence  representation. 
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r  <  X  r  <  X  r  <  X 

a.  Time  scale  effects 


r  >  X  r  <  X 

b.  Length  scale  effects 

Figure  20.  Particle  interaction  with  a  turbulent  eddy  -  effects 
of  time  and  length  scales. 
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Figure  24.  Experimental  and  predicted  particle  dispersion  for 
the  experiment  of  Snyder  and  Lumley  (Ref.  25). 
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Figure  25.  Experimenlal  and  predicted  particle  turbulent 
velocity  decay  for  the  experiment  of  Snyder  and 
Lumley  (Ref.  25). 
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b.  x/M  =  124 

Figure  26.  Predicted  particle  response  to  grid-generated  turbulence. 
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Figure  27.  Product  of  tbe  Stokes  number  and  the  turbulence 
time  scale  for  grid-generated  turbulence  for  x/M  of 
80  and  124. 
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Figure  28.  Flow  structure  for  a  simple  subsonic  jet. 
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a.  Mixing  region 


b.  Main  region 

Figure  31.  Radial  distribution  of  radial  velocity. 
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b.  Main  region 

Figure  32.  Radial  distribution  of  turbulent  stresses. 
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Figure  35.  Measured  and  predicted  particle  velocities  on  the 
jet  centerline  for  20-^m  fly-ash  particles. 

uo  Calc.  Exper. 


Figure  36.  Measured  and  predicted  particle  diffusivity  on  the 
jet  centerline  for  20-jun  fly-ash  particles 
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Figure  39.  Product  of  the  Stokes  number  and  the  turbulence 
time  scale  for  a  subsonic  axisymmetric  jet  for  x/D 
of  8.57  and  12.0. 
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b.  x/D  =  8.57,  r/D  =  0.34 

Figure  41.  Approximations  and  theoretical  particle  response 
for  a  subsonic  jet. 
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Table  1.  Particle  Characteristics  for  the  Experiment  of  Snyder  and  Lumley  (Ref.  25) 


Particle 

Diameter,  m 

Density,  g/cc 

Stokes  No. 

Hollow  Glass 

46.5  ±  3.5 

0.26 

579.0 

Solid  Glass 

87.0  ±  5.0 

2.50 

17.24 

Corn  Pollen 

87.0  ±  5.0 

1.00 

43.08 

Copper 

46.5  ±  3.5 

8.90 

16.95 
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NOMENCLATURE 

A  Coefficient  in  Eq.  (34) 

A0  Particle  acceleration  number 

Ap  Particle  surface  area 

S  Coefficient  in  Eq.  (34) 

Cd  Particle  drag  coefficient 

Cdo  Stokes  drag  coefficient 

Cjj  Coefficients  defined  in  Eq.  (67) 

D  Jet  diameter 

d  Particle  diameter 

Fi  Inertia  correction  for  Stokes  drag 
F2  Compressibility  correction  for  Stokes  drag 

F3  Rarefaction  correction  for  Stokes  drag 

F4  Shape  correction  for  Stokes  drag 

Fl  Particle  lift  force 

fi  Defined  by  Eqs.  (15),  (19),  (23) 

(2  Defined  by  Eqs.  (16),  (20),  (24) 

O  Defined  by  Eq.  (40) 

g  Gravitational  acceleration 

h  Defined  by  Eq.  (41) 
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J  Jacobian 

K  Defined  by  Eq.  (9) 

Kn  Particle  Knudsen  number 

La  Turbulence  integral  length  scale 

Im  Defined  by  Eq.  (93) 

l  Prandtl’s  mixing  length 

Mr  Particle  relative  Mach  number 
m  Grid  mesh  spacing 

mf  Displaced  fluid  particle  mass 

mp  Particle  mass 

Nf  Stokes  frequency  number 

Ns  Stokes  number 

Re2  Particle  relative  Reynolds  number  evaluated  with  conditions  behind  a  normal  shock 

ReR  Particle  relative  Reynolds  number 

Rzi  Random  variate 

r  Radial  direction 

f  r/D 

S  Density  ratio,  Gp/gf 

T  Temperature 

TEi  Turbulence  integral  time  scale 
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t 

Time 

u 

Fluid  velocity 

U 

E  ,  / 

i  Uj  Uj 

V 

Particle  velocity 

Particle  volume 

w 

Particle  relative  velocity 

X 

Axial  direction 

Xi 

Spatial  coordinate 

X 

x/D 

Xm 

Distance  normal  to  oblique  shock 

Y(t) 

Particle  position  at  time  t 

Of 

Defined  by  Eq.  (10) 

0 

Defined  by  Eq.  (1 1) 

y 

Ratio  of  specific  heats 

Ah 

Basset  term  correction  coefficient 

E 

Turbulent  dissipation 

Cpr 

Particle  diffusivity 

V 

Kolmogrov  length  scale 

*n 

(?-0.5)/x 

Vl 

r/x 

e 

Defined  by  Eq.  (53) 
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X 

Taylor  microscale 

p 

Absolute  viscosity 

V 

Kinematic  viscosity 

Cii 

Defined  by  Eq.  (14) 

Qf 

Fluid  density 

6p 

Particle  density 

T 

Time 

<t> 

Defined  by  Eq.  (14) 

Particle  shape  parameter 

a 

Defined  by  Eq.  (8) 

Subscripts 

i 

Spatial  location 

c 

Fluid  particle  location 

f 

Fluid 

P 

Particle 

X 

Particle  position 

0 

Reference  condition 

Superscript 

i 

Iteration  number 
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